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An excess in 7-rays emanating from the galactic centre has recently been observed in the Fermi- 
LAT data. We investigate the new exciting possibility of htting the signal spectrum by dark matter 
annihilating dominantly to a Higgs-pseudoscalar pair. We show that the ht to the 7-ray excess for 
the Higgs-pseudoscalar channel can be just as good as for annihilation into bottom-quark pairs. 
This channel arises naturally in a full model such as the next-to-minimal supersymmetric Standard 
Model (NMSSM) and we hnd regions where dark matter relic density, the 7-ray signal and other 
experimental constraints, can all be satished simultaneously. Annihilation into scalar pairs allows for 
the possibility of detecting the Higgs or pseudoscalar decay into two photons, providing a smoking- 
gun signal of the model. 


I. INTRODUCTION 

Cosmic-ray experiments are a promising way to search for dark matter (DM). In particular, the satellite-based 
experiment Fermi-LAT is able to measure the 7 -ray sky with unprecedented precision. If DM annihilates into 
photons with energies from 20 MeV to 300 GeV, an imprint can be left on these measurements. Intriguingly, an 
excess from the galactic centre (GCE) consistent with the range of DM density profiles indicated by observations 
and simulations of structure formation, has been identified in the Fermi-LAT data [1-3]. 

Taking the estimated uncertainty in the high-energy tail of the spectrum into account, DM annihilating to 
Higgs pairs close to threshold provides a good fit to the GCE [4, 5]. In this work we analysed the new possibility 
of fitting the signal via DM annihilation into a Higgs and a pseudoscalar and study well-motivated DM models 
in which these channels arise naturally. We find that an acceptable fit is obtained for pseudoscalar masses up 
to around 150 GeV. In particular, in the region where the pseudoscalar is lighter than the Higgs, the fit is 
improved compared to the annihilation channel into Higgs-pairs [ 6 ]. Eor a sufficiently light pseudoscalar, the fit 
becomes even better than that for the 66 -channel. 

Within the context of supersymmetric models, the light pseudoscalar required for the Higgs-pseudoscalar 
channel is difficult to obtain in the minimal supersymmetric standard model (MSSM) due to collider and flavour 
constraints [7, 8 ], [9]. Given these problems in the MSSM, we consider the NMSSM, where a light pseudoscalar 
is easier to obtain. Eurthermore, additional contributions to the Higgs quartic coupling in the NMSSM alleviate 
the need for large stop-sector soft masses to increase the Higgs mass, improving the naturalness of the model [10]. 
Both the Higgs and the pseudoscalar can decay into two photons and although we And that the 7 -ray line from 
pseudoscalar decays is distinguishable from the continuum only in an optimal-case scenario, more sensitive 7 -ray 
experiments in the future may be able to detect it. Searching for these lines in the 7 -ray spectrum could provide 
a smoking-gun signal of the new channels. 


II. THE 7-RAY excess FROM THE GALACTIC CENTER 


We next discuss the annihilation channels to a Higgs-pseudoscalar pair {ha) in the context of a two-Higgs- 
doublet model of type II with an arbitrary number of additional singlets. This includes the NMSSM which we 
focus on later. We fix the Higgs mass to 125 GeV and its couplings to SM values, as implemented in PYTHIA 
8.201 [11]. We denote the pseudoscalar as a and set tan/3 = 3 to fix its couplings. Furthermore, we assume 
that neither a nor h can decay into other scalars. 

We have performed our own fits to the reduced spectra of [3]. For completeness, we have also performed fits 
for the 66 - and /i/i-channels. We assume a generalised NEW profile, 


p{r) = PQ 



f 1 + '^q/Rs 
\ 1 r / Rs 


3-7 


( 1 ) 


with slope parameter 7 = 1.26, scale radius Rs = 20 kpc and the DM density pQ = 0.4GeV/cm^ at the 
radial distance of the sun from the galactic centre Tq. We use the prompt photon spectrum, dN^jdE, for 
annihilations into 66 from PPPC4MID [12] including electroweak corrections [13] and And good agreement with 
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FIG. 1: (a) Spectrum of the GCE from [3] and best-fit spectra for DM annihilation to ha for different values 
of rUa. (b) Best-fit regions in the DM mass and cross section for rria = 63, 90,120 and 150 GeV, with contours 

delimiting the la-, 2a- and 3a-regions. 


our own simulation using PYTHIA 8.201 [11]. We only consider the dominant decay channels of a to bb, r+r“, 
cc, photons and gluons, and simulate the resulting prompt photon spectra for hh, ha final states using PYTHIA 
8.201 [11]. Note that these spectra are unaffected by a possible singlet admixture of the Higgs or pseudoscalar. 
Indeed, such an admixture reduces their total decay widths, but the branching fractions remain unaffected (to 
leading order). The differential fiux measured by Fermi-LAT is given by. 
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with the line-of-sight integral fi^^ds over the squared DM density. The coordinate r is centred on the galactic 
centre and can be expressed as r^{s,'ip) = Tq + 5 ^ — 2rQSCos'ip, where s is the line-of-sight distance and is 
the aperture angle between the axis connecting the earth with the galactic centre and the line-of-sight. If DM 
annihilates into multiple final states, the different fluxes are summed over. 

We use the reduced spectrum of the GCE from Ref. [3] and the corresponding covariance matrix of the flux 
uncertainties including statistical and systematic errors which is publicly available. Eor simplicity we keep the 
pseudoscalar mass rua fixed and perform a two-parameter fit in the DM mass and annihilation cross section. 
In Eigs. 1 we show the resulting best-fit spectra from DM annihilation together with the spectrum of the GCE 
for different values of and the different annihilation channels. The smallest in these figures is chosen 
such that Higgs decays to pseudoscalars is kinematically forbidden, i.e. rUa > rrihl^. The salmon-colored boxes 
depict the empirical model systematics [3], the error bars correspond to the statistical errors, and the yellow 
boxes are the combination of the statistical errors, empirical model systematics and other systematics. 

Notice that, similar to the hh-channel, the spectra for the ha-channel in Eig. 1 have a peak at energies m/j,/2 
which is produced from on-shell decays of the Higgs to two photons. The peak is less pronounced than for the 
hh-channel because there is only one Higgs in the final state. There is only mild line broadening for the best-fit 
masses, since the hh- or ha-pair is produced close to threshold. Notice also that the spectra have no visible 
peak at energies ma/2 from pseudoscalar decays to two photons. Indeed, the peak in Eig. 1 for the ha-channel 
is purely due to Higgs decays. Using the Eermi-LAT limits on y-ray lines [14], it was estimated that the line 
strength from the decay hh ^ dy is just below current limits and may be detected (or excluded) in the near 
future [3]. In contrast, the intensity of the line produced by pseudoscalar decays is too weak to be in tension 
with line searches. We find that only in the best-case scenario of rria ^ 150 GeV and tan^d ^ 1 is the peak 
barely distinguishable from the continuum. More sensitive y-ray experiments may be able to detect the photon 
peak for lower rUa and larger tan [3 in the future. 

In Eig. 1 we also show the best-fit regions in the DM mass and cross section for the ha-channel. Eor 
comparison, we also show the best-fit regions that we find for the bb- and hh-channels. Similar to [4, 5], we 
assume a multiplicative astrophysical-uncertainty factor A for our best-fit cross-sections. This factor takes into 
account the uncertainties in the local DM density p©, the scale radius Rs and the slope parameter y. We have 
used the same reference values for pQ, Rg and y for our fits as [4] and we therefore use their estimate for the 
range of the astrophysical-uncertainty factor, A G [0.17,5.3]. Notice from Eig. 1 that the best-fit regions for 
the ha-channel lie very close to threshold. The best-fit cross sections are fairly independent of the pseudoscalar 
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channel 

ma [GeV] 

"Idm [GeV] {av] 

>0 [10 ^®cm^/s] 

Xmin 

p-value 

66 


49.61®:' 

1 r:+0.3 
-'-•^-0.2 

24.5 

0.32 

hh 


125.012:^0 

A 2 + 0-® 

30.0 

0.12 


63 

94.0112 

o ^+0.7 

22.4 

0.43 

ha 

90 

107.5111 

q 7+0.7 
'J- ‘ -0.7 

25.3 

0.28 


120 

122.5111 

3.811? 

30.3 

0.11 


150 

137.511^ 

4.011? 

36.0 

0.03 


TABLE I: 


Best-fit values and Icr-regions from our fits for the bb, hh, ha annihilation channels. 
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FIG. 2: Feynman diagrams contributing to neutralino annihilation to ha-pairs. Note that there is a i^-channel 

diagram in addition to the t-channel one. 


and DM mass. In Table I, we show the values and la-regions for the best-fit DM mass and cross section and 
the associated p-values for the different annihilation channels. The smaller rria lead to a better fit, i.e. a 

smaller and thus a larger p-value. In particular, for rua ^ 120 GeV, the fit for the ha-channel is better than 
for hh final states and for rUa = 63 GeV the fit becomes better than for bb final states. It can also be seen that 
the best-fit cross sections for the scalar channels are larger than for the 66-channel. 

Let us briefly mention that there is no dedicated analysis for annihilation into hh and ha in the preliminary 
dwarf-spheroidal galaxy searches. However, 6-quarks are still the dominant decay product for these channels. 
We obtain the rough estimate (an)o ^ 3 x 10“^^cm^/s for DM with ^ 100 GeV annihilating into hh or 
ha close to threshold, depending relatively weakly on rria in the range of interest here. 


III. ANNIHILATION INTO ha IN THE NMSSM 

The annihilation of neutralinos is p-wave suppressed if the final state is even under CP, such as for hh or 
aa [15]. Therefore it is only through the ha-channel that neutralino annihilation into scalars can account for 
the GGE. The graphs that contribute to the corresponding cross section are the s-channel exchange of the two 
pseudoscalars and the t/rt-channel exchange of the neutralinos. These are shown in Fig. 2. We have performed 
two random scans using NMSSMTools 4.4.0 [16] for the ha-channel, one optimised for bino-like LSPs and the 
other for singlino-like LSPs. The masses of the squarks, sleptons, gluino and wino are fixed at 2 TeV and the 
remaining free parameters are scanned over the ranges shown in Table 11. Gollider constraints and direct and 
indirect detection constraints are all satisfied. Here A\ is partly determined by the requirement that the singlet 
admixture to the Higgs is small. In Fig. 3, we show a scatter plot of the annihilation cross-sections during 
freeze-out and at late times. The approximate degeneracy of these cross sections reflects the dominance of the 
XX ^ ha process and the absence of a large special enhancement such as a resonance. In Fig. 3, we also show 
a scatter plot of the cross section for spin-independent DM-nucleon scattering versus the DM mass. The LUX 
collaboration will test a considerable portion of the singlino-like points and almost the entirety of the bino-like 
points from our scan of the ha-channel. The vast majority of points with a singlino-like LSP would then be 
probed by XENONIT, although we find a few such points that would evade even the projected LZ limits. 


[GeV] [GeV] 

Ateff [GeV] Ml [GeV] 

A n 

tan 13 

bino-like [-50,50] [-100,100] 

singlino-like [-50,50] [-100,100] 

[-300,-100] 

[-600,-200] 

[60,170] 

2000 

[0.6,1.4] [0.1,1.6] 
[0.6,1.4] [0.05,0.5] 

[2,5] 

[2,5] 


TABLE H: Parameter ranges for the two random scans for the ha-channel. The first and second line are for 
the scan optimised for bino-like and singlino-like LSPs, respectively. 
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FIG. 3: The left panel (a) shows the annihilation cross section at the present time, {(Jha'^)oi and at freeze-out, 
The right panel (b) shows the spin-independent cross section for scattering of LSPs off nucleons 

versus the LSP mass. 


IV. CONCLUSIONS 

We have investigated the possibility of explaining the 7 -ray excess at the center of the galaxy via DM 
annihilation into a Higgs and a pseudoscalar, and showed that a good fit of the photon spectrum can be 
accomplished. We have also shown that this possibility arises naturally in the context of the NMSSM where 
regions exist that satisfy all experimental constraints. 
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